c-axis conductivity in the normal state of cuprate superconductors 
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The c-axis optical conductivity and d.c. resistivity are calculated within the t-J model assuming 
that the interlayer hopping is incoherent. Use is made of numerical results for spectral functions 
recently obtained with the finite-temperature Lanczos method for finite two-dimensional systems. In 
the optimally doped regime we find an anomalous relaxation rate r" 1 oc u) + £ c T and p c (T) oc p a b(T), 
suggesting a common relaxation mechanism for intra- and interlayer transport. At low doping a 
pseudo-gap opening in the density of states appears to be responsible for a semimetallic-like behavior 
of P c(T). 

PACS numbers: 71.27.+a, 72.15.-v 
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One of the most striking characteristics of cuprate su- 
perconductors is the anisotropy of their structures. The 
Cu02 planes common to all high-T c materials clearly de- 
termine most of the normal-state electronic properties 
which are quite anomalous but at the same time rather 
universal within the whole family of cuprates. The quan- 
tity which most evidently displays the anisotropy of a 
particular material is the optical conductivity a(u>,T) 
and the corresponding d.c. resistivity p{T) [fil[|. The 
ratio C(T) = p c j Pah of the c-axis resistivity p c to the 
in-plane resistivity p a b for temperatures T just above T c 
in optimum-doped materials ranges from £ ~ 100 for 
YBa 2 Cu 3 7 _ 5 (YBCO), ( ~ 1000 for La 2 _ a Sr !r Cu04 
(LSCO) to C > 10 5 in Bi 2 Sr 2 CaCu 2 8 . 

In spite of such large quantitative differences there are 
common qualitative features even in the c-axis trans- 
port, in particular regarding the variation with doping, 
(a) Optimum-doped YBCO and overdoped LSCO show 
metallic-like p c {T) ~ a+bT for T > T c . At the same time 
the optical conductivity cr c (u>) shows a fall-off with oj al- 
though far from the usual Fermi-liquid (FL) Drude-type 
behavior, (b) For lower doping the variation is gener- 
ally semimetallic-like with dp c /dT < in the interval 
T c < T < T* . The crossover temperature T* showing up 
e.g. as a kink in the in-plane p a b(T), in the Hall constant 
Rh(T) etc. ||, decreases with doping and merges with 
T c at the optimum doping. Within the low doping regime 
a c (uj) is nearly flat in a wide frequency regime for higher 
T while a pseudo gap (PG) starts to emerge for ui < lo* 
for lower T < T pg < T*. 

The c-axis conductivity has attracted the attention of 
theoreticians since the first measurements on cuprates 
W . In spite of open challenging questions of anomalous 
planar properties, there seems to be an agreement in the 
conclusion that the c-axis conductivity generally cannot 
be explained as a coherent interlayer transport (jj^J . The 
main indication is that the d.c. a c is mostly well below 
the Mott minimum metallic conductivity ||, with the 
possible exception of optimum-doped YBCO and over- 
doped LSCO fy. In other words, estimates of the mean 



free path in the c-direction appear much shorter than 
the interlayer distance cq, i.e. l c {T) <C cq, when one as- 
sumes the usual Boltzmann approach for metals. The 
absence of coherent c-axis conduction is a consequence 
of weak interlayer hopping matrix element t c but also 
of a strong intralayer scattering. Based on the concept 
of dynamical detuning, Leggett |0j thus derives the con- 
dition for the incoherent conduction t c <C l/r Q ~ k-&T 
taking the relaxation-rate according to anomalous pla- 
nar p a b(T) and a a b(uj). The incoherent conductivity has 
been in more detail studied on the example of coupled 
fermion chains p] representing Luttinger liquids. Still 
there have been so far no real attempts to derive more 
explicitly the c-axis conductivity for models relevant to 
cuprates. 

In this paper we present the theoretical analysis of 
a c {ijj) within the simplest microscopic model incorporat- 
ing both strong electron correlations in the CuO planes 
leading to antiferromagnetism (AFM) in undoped mate- 
rials, and a weak interlayer coupling. Within each layer 
we consider the planar t-J model 
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+R.c.)+J^2(Su-Si j -jni i ni j ), (1) 
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where i,j refer to planar sites on a square lattice within 
the l-th layer and cu s ,cj is represent projected fermion 
operators forbidding double occupation of sites. As ap- 
propriate for cuprates we assume J — 0.3 t. We allow for 
the hopping between layers, 



Q+li 



H.c). 



(2) 
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The planar t-J model is still a challenge for theoreti- 
cians, nevertheless there is a growing consensus that it 
contains essential ingredients for an explanation of the 
anomalous normal-state properties of cuprates, domi- 
nated by effects of strong correlations. The lack of reli- 
able analytical methods in the latter regime led to nu- 
merous numerical studies M). Using the novel finite- 
temperature Lanczos method (FTLM) M it has been 
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shown that several thermodynamic quantities as well re- 
sponse functions match well as the experimental data on 
cuprates. Most relevant for the present study are the re- 
sults on the planar conductivity cr a b(u>) an d the single- 
particle spectral functions (SF) A(h,u>) JlCfl. 

In the limit of weak interlayer hopping t c <C t it is 
straightforward to evaluate the dynamical c-axis conduc- 
tivity <j c (uj) using the linear response theory. Express- 
ing the c-axis current correlation function (j c (t)j c ), with 
jc = e ot c coY,u s ( ic iis c Ws +H.c), in terms of the planar 
SF A(k, u>), we arrive at 



a c (u) = ^ I dw'[f{J)-f{uj' + w)] 

UJ 



x J2 A(k, (j')A(k, uj' + uj), 

k 



N 



(3) 



where cr° = e\t\c^lTia^t 2 is a characteristic c-axis con- 
ductivity scale, N is the number of sites and f(ui) 
' 1 + exp(w/fcBT)] _1 . The approximation assumes the 
independent electron propagation in each layer, and is 
justified for t c <C i. Note that the interlayer-hopping 
term, Eq. (Q), conserves the QP momentum k, and in 
this respect the treatment is analogous to the one of 
the transverse transport in weakly coupled Luttinger 
chains || . This could be compared with a related prob- 
lem of interlayer hopping with random matrix elements 
(t c )u, where only the planar density of states (DOS) 
j\f(u>) — 2/N uj — fi) enters, with /i denoting the 

chemical potential. The corresponding expression for the 
\i-nonconserving (J™(oj) is obtained by the replacement 



xirt 2 Af(i.i + oj')j\f(fi + oj' + oj), 



(4) 



appearing e.g. in disordered systems p], together with 
a substitution for t c in the definition of er|? with some 
average t c . In cuprates both alternatives have a justifica- 
tion, since the disorder introduced by dopands residing 
between layers modifies the hopping elements. Hence one 
can expect that the actual conductivity is a linear com- 
bination of <7 C and cr™. 

The knowledge of planar A(k, uj) and Af(fj, + uj) should 
thus suffice for the evaluation of <r c (uj) . In cuprates, how- 
ever, SF are available with sufficient resolution only in 
the hole part uj < (in principle) via the angle-resolved 
photoemission (PES) and the DOS via the angle inte- 
grated PES @. Within the t-J model (and the closely 
related Hubbard model) A(k, uj) and jV(uj) have been 
studied mostly numerically 0, applying the exact diag- 
onalization (ED) of small systems ]lo| and the quantum 
Monte Carlo method Jl4| . Results reveal at intermediate 
doping a large Fermi surface (consistent with the Lut- 
tinger theorem) and a quasiparticle (QP) dispersion ek 
similar to but reduced in bandwidth relative to the free 



tight-binding electrons. Only recently, via the FTLM 
a reliable evaluation of the corresponding self en- 
ergy S(k, uj) was possible, thus allowing for a study of 
low-T QP properties. The SF and a a b(uj) both reveal in 
this regime the anomalous behavior consistent with the 
marginal Fermi liquid (MFL) concept p5| , i.e. the effec- 
tive transport relaxation rate as well as the QP damping 
appear to follow 1/T a b(uj) oc S"(k, oj) oc \uj\ + £, a bT. 

In this paper we use numerical results for A(k, oj) and 
jV(uj), as obtained using FTLM on systems with N — 
16, 18, 20 sites for various hole dopings Ch = Nh/N. SF 
are available for the whole range of Nh for N = 16, Nh < 
2 for TV = 18, and N h < 1 for N = 20 0J|. One should 
keep in mind the restriction that FTLM results become 
dominated by finite-size effects for T < Tj s . Tj s is clearly 
size and doping dependent. The lowest T can be reached 
in the intermediate regime Ch ~ 0.2 where Tf s ~ 0.1 t, 
while Tf s increases both for lower and higher doping. To 
make contact with experiments, we note that t ~ 0.4 eV, 
i.e. our numerical results apply to T > 450 K, while for 
lower T they can give some qualitative indications. 

To set the frame for further discussion let us first 
present in Fig. 1 the DOS W(m + w) for the two low- 
est nonzero doping levels Ch = 1/20 and Ch = 2/18 for 
several T < J. The DOS at larger doping ct ^ 0.15 
(not shown here) E3] is T-independent and featureless 
at oj ~ 0, as expected for a Fermi liquid away from van 
Hove singularities. At low doping the DOS in Fig. 1 is 
also structurelesss for T > J, but here this reflects the 
incoherent character of QP due to large damping. Upon 
lowering T a gradual transfer of weight in the underdoped 
samples from above fi to oj < is observed, leading to the 
formation of a PG at uj ~ 0. The deepening is more pro- 
nounced for lowest dopings, e.g. = 1/20 in Fig. 1(a), 
where the PG energy scale is well below t and appears to 
be A pg ~ J, reflecting the onset of the short-range AFM 
order. 
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FIG. 1. The DOS A% + uj) at various T < J for hole 
concentrations: (a) cu = 1/20 and (b) Ch = 2/18. 
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The above development of the DOS as well as the be- 
havior of SF and corresponding (anomalous) self energies 
should determine the behavior of optical conductivity 
a c (uj). In Fig. 2 we show a c (oj) for two dopings Ch = 3/16 
and c h = 1/20 for T < J, as calculated from Eq. (g). At 
intermediate doping Ch = 3/16 a c (oj) in Fig. 2 exhibits a 
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FIG. 2. a c (uj) at several T for dopings: 
c h = 1/20. 



Ch = 3/16 and 



central peak, sharpening at lower T. This can be ex- 
plained with the properties of A(k, oj) in this regime, con- 
sistent with a large Fermi surface and a QP dispersion ek 
corresponding to a tight-binding band. The QP damp- 
ing is however large and MFL-like S"(k, oj) cx \oj\ + £ Q t,T, 
but only weakly k-dependent. Inserting the latter into 
Eq. (||) it is easy to see that the resulting a c {oj) also fol- 
lows the MFL-behavior, i.e. one can represent it in the 
generalized Drude form with an effective relaxation rate 
1/t c oc |o;| + £ C T '. The vanishing rate for oj,T — > is 
clearly the consequence of the k-conservation in the in- 
terlayer hopping. On the other hand, Eq. (^) would yield 
a qualitatively different result, i.e. an almost flat and T- 
independent cr™(u;) due to the structureless DOS j\f(oj) 
for oj ~ 0. 

The corresponding results for low-doping = 1 /20 are 
also shown in Fig. 2. In contrast, we find rather weakly uj- 
and T-dependent a c (u>), both for the k-conserved and the 
k-nonconserved approximation. Quantitative agreement 
of both approaches indicates that in this regime we are 
dealing with a strongly reduced effective dispersion ek oc 
J but still large QP damping, hence weakly k-dependent 
A(k,uj) ~ N(p + oj), at least for T > Tf s . It seems 
puzzling how rather constant cr™(w) in Fig. 2 for Ch — 
1/20 can be compatible with a well pronounced PG at 
oj ~ in the DOS (see Fig. 1) at low T < T* . These 
facts are however reconciled by observing that the Fermi 
functions in Eqs. (g^) have an effective width ~ 4T, 
hence down to reachable Tf s they smear out both the 
QP dispersion as well as the A pg . Nevertheless, there is 
some signature of a PG opening on a scale oj ~ J in a c (oj) 



for lowest T, in qualitative agreement with experimental 
data |16|. 

In Fig. 3 we present the d.c. resistivity p c (T) = 
l/a c (oj = 0, T) for a wider range of Ch- In the regime 
of intermediate and even higher doping Ch > 3/16 p c {T) 
is metallic-like for all T. In particular from the previous 
arguments relating the c-axis transport to the MFL-like 
planar relaxation one expects that pc(T) oc T for T < J. 
From Fig. 3 we realize that the latter behavior is re- 
ally restricted to T < 0.2 t. It should be stressed again 
that in this regime weak interlayer disorder as implied by 
Eq. (|3|) is essential for the resulting linear-in-T behavior 
of p c (T), whereas p c (T) in the case of k-nonconserving 
interlayer hopping (not shown in Fig. 3) yields an almost 
T-independent p c {T). 




FIG. 3. pc(T) vs. T for various Cfe. Solid and da- 
shed lines represent results within the k-conserving and 
k-nonconserving approximation, respectively. 

For the lowest c h = 2/18 and Ch — 1/20 p c and 
yield the same qualitative behavior, p™(T) being system- 
atically higher. While p c (T) is increasing for T > J, t as 
generally expected for incoherent transport, the relevant 
question is its behavior for T < J. In Fig. 3 we notice 
that in this regime dp c /dT ~ 0, signalling the onset of a 
nonmctallic-like behavior in agreement with experimen- 
tal results for cuprates in the underdoped region Jl7t ■ 

When we comment on the relation of our results with 
measurements on cuprates, the most natural application 
is to LSCO having the simplest layered structure among 
cuprates. The presence of several nonequidistant layers 
per unit cell, and possibly chains as well, in other ma- 
terials complicate the interpretation of the c-axis trans- 
port considerably. In Fig. 4 we present the resistivity 
ratio C(^) = Pel Pab as a function of T for different dop- 
ing levels, where results for p a b are taken from previous 
FTLM studies |J . In the saturation regime T > J we set 
C ~ 150 as obtained from experimental data for LSCO 
e.g. by Nakamura and Uchida JTt]] and Kao et al. Ej| 
from which we can estimate t c . Taking the standard 
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value for co/ao — 3.5 in LSCO Q] we get t c ~ 0.03i, 
or i c ~ 12 meV. 

A distinctive feature is that £(T) is almost indepen- 
dent of c/j and T, but for the lowest T < J presented. 
This indicates that planar resistivity and c-axis resistivity 
are related, at least for T > Tf s , pointing to a common 
current relaxation mechanism in all directions. In fact, 
if we calculate within the same decoupling scheme the 
planar conductivity er^ c (w) and the related p^(T), the 
actual a a b{w) is reproduced remarkably well, as can be 
concluded from ( ab (T) = p*f / p ab in Fig. 4. For T - J 
we detect the onset of a different behavior, although the 
ratio is less affected than p a b(T) and p c {T) separately. 
It is however indicative that the increase of C(T) is most 
pronounced for lowest c/j. Our results for ((T) are in 
a remarkable qualitative agreement with experimental 
data p). 



spin correlations, persisting up to the optimum doping. 
It is also evident from our results that this does not lead 
directly to a well pronounced PG in <j c {uj) at the same T. 
Eventually this can happen only for T < T pg < T* , which 
in our systems is hardly reachable due to Tf s ~ 0.2 1 
at low doping. One simple explanation can be given in 
terms of the Fermi- function broadening ~ 4T. Another 
source of suppression of the PG in calculated <j c {oj) is the 
shift of the PES leading edge towards p (see Fig. 1), thus 
compensating in part for the deepening of the PG in the 
DOS at lowest T. This results in a much weaker dop- 
ing dependence of p c down to T ~ Tf s , not inconsistent 
with experimental data []l7f which show semiconductor- 
like upturn only at lower T . Still it remains the subject 
of future studies to clarify whether the small PG scale in 
underdoped cuprates is directly related to the T* scale 
or is of different origin. 
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FIG. 4. pc/ Pab vs. T for various Ch (left scale) and 
Cab(T) = ptr/pab (right scale). 

In conclusion, let us briefly summarize our main re- 
sults on the c-axis conductivity. In our analysis a c (u>,T) 
is obtained from (numerically) known SF and DOS for 
the planar t-J model supplemented with a small (doping 
independent) interlayer hopping t c leading to the incoher- 
ent c-axis transport. A qualitative agreement with exper- 
iments in LSCO, e.g. for the anisotropy C(T), gives sup- 
port to such a minimum model. Our results also confirm 
the general experimental observation that p c (T) changes 
from a metallic one to a semimetallic one dp c (T) / dT < 
by decreasing doping. It should be however noted, that 
in the overdoped regime LSCO shows already indications 
of an enhanced coherent-like c-axis transport which is be- 
yond our analysis assuming t c <C t. 

Let us finally point on open questions, mainly related 
to the onset of a PG in ct c (u)) and even semiconductor- like 
dp c (T)/dT <C as found experimentally in underdoped 
cuprates In our analysis the indication for such a 
development comes from a PG in the DOS, which at low 
doping starts to emerge for T < T* < J, and is com- 
patible with recent PES experiments [Q. Clearly, this 
phenomenon is related to the onset of short-range AFM 
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